Rheophysics of dense granular materials : 
Discrete simulation of plane shear flows 
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We study the steady plane shear flow of a dense assembly of frictional, inelastic disks using 
discrete simulation and prescribing the pressure and the shear rate. We show that, in the limit of 
rigid grains, the shear state is determined by a single dimensionless number, called inertial number 
J, which describes the ratio of inertial to pressure forces. Small values of / correspond to the quasi- 
static regime of soil mechanics, while large values of I correspond to the collisional regime of the 
kinetic theory. Those shear states are homogeneous, and become intermittent in the quasi-static 
regime. When / increases in the intermediate regime, we measure an approximately linear decrease 
of the solid fraction from the maximum packing value, and an approximately linear increase of the 
effective friction coefficient from the static internal friction value. From those dilatancy and friction 
laws, we deduce the constitutive law for dense granular flows, with a plastic Coulomb term and a 
viscous Bagnold term. We also show that the relative velocity fluctuations follow a scaling law as 
a function of /. The mechanical characteristics of the grains (restitution, friction and elasticity) 
have a very small influence in this intermediate regime. Then, we explain how the friction law 
is related to the angular distribution of contact forces, and why the local frictional forces have a 
small contribution to the macroscopic friction. At the end, as an example of heterogeneous stress 
distribution, we describe the shear localization when gravity is added. 

PACS numbers: 45.70.Mg, 81.05.Rm, 83.10-y, 83.80.Fg 



I. INTRODUCTION 

Due to their importance in geophysics (propagation of 
avalanches, migration of dunes, fault sliding) and in vari- 
ous industrial processes (handling of powders, granulates 
in civil and chemical engineering, food, pharmacy, tri- 
bology...), flows of granular materials are the focus of a 
large number of research, at the frontier between physics 
and mechanics 0, B IE H • I n order to predict propaga- 
tion, flow-rate, jamming..., one of the main objectives of 
these "rheophysical" studies is to determine the rhcologi- 
cal laws of those materials, based on their physical origin 
at the scale of the grains and of their interactions. One 
thus tries to express the stress tensor (and especially the 
pressure P and the shear stress S, positively counted) as 
a function of the shear rate 7 and other variables such as 
solid fraction v. 

Granular materials are extremely various depending 
on the geometry of the grains (shape and size) and the 
nature of their interactions. We shall restrict our atten- 
tion in the following to assemblies of cohesionless grains, 
slightly polydisperse, without interstitial fluid. This cor- 
responds to macroscopic grains (diameter larger than 
hundred microns) in a fluid of low viscosity like air. The 
rheology is then only dictated by transfer of momentum 
and dissipation of energy taking place in direct contacts 
between grains and with the walls. 

Depending on the conditions, these materials reveal 
various mechanical behavior, similar to elastoplastic 



solids in the quasi-static regime, to dense gazes in the 
case of strong agitation, or to viscoplastic fluids when a 
flow is provoked. This paper is devoted to this interme- 
diate regime, which is still not well understood. We shall 
first briefly recall the essential results for the two extreme 
regimes, quasi-static and collisional. 



A. Quasi-static regime 



'Electronic address: Corresponding~author : chevoir@lcpc.fr 



Dense, confined granular assemblies in extremely slow 
shear flow (7 = 10~ 4 s -1 for sand) are usually de- 
scribed as solids abiding by elastoplastic constitutive 
laws [EIESI3) or more general, incremental ones P.ITo|. 
A typical experiment consists in imposing a very slowly 
growing deviator (shear) strain to a sample prepared in 
mechanical equilibrium under an isotropic pressure. The 
material response is then independent of physical time, 
and involves volumetric strain, i.e., dilatancy : loose 
samples contract, whereas dense ones expand. In dense 
samples, the shear deviator stress reaches a maximum 
(around shear strains of order 5%) and then decreases to 
a plateau. In a loose system, it increases monotonically to 
the same plateau value under the same confining stress. 
Meanwhile, the solid fraction approaches a limit which 
does not depend on its initial value either. It is therefore 
generally accepted that after a large enough shear strain 
(of order 10%), the material has reached a certain attrac- 
tor state, called the critical state [TH , which does not 
depend on its initial arrangement. If larger shear strains 
are imposed, the critical state, with its internal structure 
adapted to shear flow in the limit of small velocities, will 
remain unaltered. The critical state is characterized by 
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an internal friction angle 0, defined as tan <j> = S/P in a 
simple shear test, and by a critical solid fraction v c . In 
the case of dense samples, the observation of the critical 
state is often precluded by the onset of shear localiza- 
tion, whereby the corresponding microstructure is lim- 
ited to the material that is continuously sheared within 
thin bands. 

Discrete numerical simulations have provided some ad- 
ditional information on the internal structure of granular 
materials as they approach their critical state. Thus the 
critical state is associated with specific values for coordi- 
nation numbers and distribution of contact orientations 
("fabric") [T^- Simulations also contributed to unravel 
the influence of contact elasticity on the macroscopic be- 
havior. In the approach to the critical state, strains are 
essentially due to rearrangements, as contact networks 
break and repair. Macroscopic strains, with typical pres- 
sure levels and contact stiffnesses, as soon as they reach 
the 10 -4 range, originate in geometry changes and not in 
contact deformability 0, . As a sample approaches 
its critical state, the gradually increasing anisotropy of 
the contact network allows it to support increasing de- 
viator stresses ^3|. Contact elasticity is essentially ir- 
relevant in that regime, which can be modelled with 
rigid grains [l^. The internal friction angle 4> an d the 
critical state solid fraction v c are found to depend only 
negligibly on the confining pressure P. The slight de- 
pendence reported in experiments on sands ^2] is likely 
due to features of the contact law, such as plasticity or 
breakage, which are usually not introduced in numerical 
simulations. In other words, the continuously sheared 
material can be modelled with the Coulomb criterion, 
S = Ptan0. If the direction of the imposed strain rate 
changes (e.g., upon reversing the flow), then the system 
will leave its critical state and gradually ap pro ach an- 
other one, with different contact orientations [l3j . 

Quasi-statically deformed granular systems are one 
limit case for the steady plane shear flows investigated 
in the present paper. Therefore, the results will be pre- 
sented in the sequel in terms of an effective friction co- 
efficient /i* = S/P, which should come close to tan0 
in the limit of slow motion. Likewise, as appropriate for 
quasi-static granular rheology, we chose to control the lat- 
eral pressure rather than keeping the solid fraction fixed. 
Such options are further commented in the following. 



B. Collisional regime 

In the dilute limit and/or for strong agitation, the 
grains interact through binary, instantaneous, uncorre- 
cted collisions. Then, the generalization of the kinetic 
theory of dense gazes to slightly inelastic grains [17], 
El El EH allows a hydrodynamical description. The 
stress components depend on the solid fraction and on 
the velocity fluctuations Sv (which square defines the so- 
called granular temperature). In the two-dimensional ge- 
ometry which we shall study in the following, the stress 



components are homogeneous to a force divided by a 
length. For an assembly of disks of diameter d and mass 
to : 

P = Fp(v)m(8v/d) 2 , 

(1) 

S = Fs(v)m(5v I 'd)j . 

Solving a flow problem requires an additional equation 
of energy in which a dissipation rate T associated to in- 
elastic collisions must be added to the usual terms : 

T = F r (v)m{6v l 'df . (2) 

The dimensionless functions are completely ex- 

pressed as functions of the pair correlation function at 
contact go(v). In the dense limit (0.2 < v < 0.67), 
Fi(v) — AiF{v), with F(v) = v 2 g n(v ). In two dimen- 
sion, g (v) = (16 - 7f)/16(l - v) 2 P3. The pre- factors 
A j ar e well-known functions of the restitution coefficient 
e 111. 

In the case of an homogeneous system, where the shear 
rate and the velocity fluctuations are uniform, the equa- 
tion of energy reduces to a balance between the work of 
the shear stre ss and t he dissipation : Sj — T. This leads 
to Sv = -jdy/ As/Ar- Consequently, the stress compo- 
nents are equal to : 

P = Gp{v)m^ 2 , 

(3) 

S = G s {v)mi 2 , 

where Gp,s{v) = Bp t sF(v) with Bp = ApAs/Ar and 
F>s = ^ Ag/Ap. We notice that the effective friction co- 
efficient is a constant (/j* = yj AsAr/A p ), and that the 
solid fraction is a function of the dimensionless quantity 
I = A /y / m/P (y — Gp 1 (l/J 2 )). This collisional descrip- 
tion is relevant in the dilute limit when the inertial effects 
dominate. 

C. Dense regime 

In the intermediate regime, the solid fraction is close 
to a maximum solid fraction, so that one speaks of dense 
flows. Then, the grains interact both through enduring 
contacts and through collisions. There exists a contact 
network more or less percolating through the material, 
which is very fluctuating in space and time |2^|. These 
flows are beyond the quasi-static regime, since the inertia 
of the grains (and so the shear rate) certainly comes into 
play. On the other hand, the assumption of binary, in- 
stantaneous, uncorrelated collisions of the kinetic theory 
is clearly in trouble. Due to the very strong correlations 
of motion and force, the theoretical description of those 
dense flows is very difficult and is still a matter of debate 
(see [24] for a recent review). 
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Advances have come in the last decade from the combi- 
nation of discrete numerical simulations and experiments 
on model materials in simple geometry, confined or free 
surface flows (annular shear cell, vertical chute, inclined 
plane, heap- flow), and in various mechanical configura- 
tions (gravity, velocity or force imposed at a wall). A 
detailed review of these works can be found in [2{| . 

Depending on the mechanical configurations, the flows 
are steady, intermittent, or even jam. A localization of 
the shear, with a width of a few grains, is also frequently 
observed near the walls or near the free surface, with ex- 
ponential velocity profiles around. However, the hetero- 
geneity of the stress distribution as well as the presence of 
walls makes the analysis of the constitutive law difficult. 



D. Organization of the article 

This is the reason why we have chosen to study this 
dense flow regime in steady homogeneous shear state. We 
have studied the simplest geometry, plane shear without 
gravity, in which the stress distribution is homogeneous 
inside the shear layer. Furthermore, we have prescribed 
both the shear rate and the pressure. Using discrete nu- 
merical simulations ("numerical experiments"), we have 
access to microscopic information, at the level of the 
grains and of the contact network, hardly measurable 
experimentally, and we are able to vary the parameters 
describing the grains and the shear state. 

Sec. [n] is devoted to the description of the simulated 
system. Then in Sec. lIIII we build a dimensionless number 
/, which, in the limit of rigid grains, describes the shear 
state and characterizes the progressive transition between 
the quasi-static and the dynamical regimes, called inter- 
mediate regime . We then show in Sec. II VI that we obtain 
homogeneous states in term of structure (solid fraction) , 
kinematics (shear rate) and stress distribution. In Sec.lVl 
we measure the evolution of two macroscopic quantities 
(solid fraction and effective friction coefficient) as a func- 
tion of I in the intermediate regime, from which we de- 
duce the constitutive law. We then study the sensitivity 
of this constitutive law on the mechanical properties of 
the grains (restitution, friction and elasticity). The next 
sections are devoted to more detailed information, first 
on the fluctuations of the grain motion (Sec. IVI|I . and 
then on the contact network (Sec. IVII|) . Then, using 
those microscopic information, we explain in Sec. IVfffl 
how the friction law is related to the angular distribu- 
tion of contact forces, and why the local frictional forces 
have a small contribution to the macroscopic friction. In 
the last Sec. IIXI we study how the previous results are 
affected when the stress distribution becomes heteroge- 
neous, in presence of gravity. For a more detailed account 
of the results, we refer to |26j . 



II. SIMULATED SYSTEM 
A. Plane shear 

The simulated system is two dimensional (Fig.^l. The 
granular material is a dense assembly of n dissipative 
disks of average diameter d and average mass m. Except 
in Sec. IIV D 21 a small polydispersity of ±20% is con- 
sidered to prevent crystallization |22|. The mechanical 
properties of those grains are described by a friction co- 
efficient fi, a restitution coefficient in binary collisions e 
and elastic stiffness coefficients k n and k t . 





FIG. 1: (Color online) Plane shear : (a) Quasi-static regime 
(I = 10~ 2 ) (b) Dynamic regime (I = 0.2,). (Black grains 
constitute the rough walls. The line widths are proportional 
to the intensity of the normal force between grains). 

The granular material is submitted to a plane shear, 
without gravity (except in Sec. IIX|) . so that the stress 
distribution is uniform. The material is sheared between 
two parallel rough walls, distant of H. One of the wall is 
fixed, while the other moves at the prescribed velocity V. 
We call x the flow direction and y the transverse direc- 
tion. Periodic boundary conditions are applied along the 
flow direction, and we call L the length of the simulation 
box (always larger than 40 grains). The wall roughness is 
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made of contiguous grains sharing the characteristics of 
the flowing grains (same polydispersity and mechanical 
properties - no rotation) . y — corresponds to the center 
of the glued grains on the fixed wall. 

An important feature of our simulation is that the nor- 
mal stress Tiyy is prescribed, as usual in practice. We 
shall see in the following (Sec. IIV B 2|) that the normal 
stress components Ti xx and ~E yy are equal, so that T, yy is 
equal to the pressure P = (T, xx + E yy )/2. The control 
of the pressure P is achieved by allowing the dilatancy 
of the shear cell along y (H is not fixed), through the 
vertical motion H of the moving wall. The evolution of 
H is given by : H = (P — P w )L/g p , where g p is a vis- 
cous damping parameter, and P w is the normal stress 
exerted by the grains on the moving wall. Steady state 
corresponds to (P w ) = P. 

The parameters of the simulated system are summa- 
rized in Tab. I (anticipating on Sec. 17771 we indicate the 
value of the dimensionless parameters g v j \/mk n and /). 



V?- — riij ■ Vij is the time derivative of the normal de- 
flection of the contact (or apparent "interpenetration" of 
undeformed disks), hij = (di + dj)/2 — \\Hij\\. Its tan- 
gential component V? = tij ■ Vij is the time derivative 
of the tangential relative displacement Sij . (Let us note 
that the definition of Sij as a scalar quantity automati- 
cally accounts for the material transport of strain in the 
contact region, as tij moves with the pair in contact). 

The normal contact force is the sum of two contri- 
butions, an elastic one N e and a viscous one N v : 



1. Viscoelasticity 

Keeping in mind that contacts have to close to transmit 
forces (Fij = if < 0) the linear (unilateral) elastic 
law reads 



n L/d H/d g p /^Jmk n I 

900 - 5000 40 - 100 20 - 100 1 6.10~ 4 - 03 

TABLE I: List of system parameters. 



Due to its interest in rheology and more specifically 
in tribology (third body) and in geophysics (sliding of 
faults at the origin of earthquakes), this plane shear ge- 
ometry has already been the subject of numerous dis- 
crete simulation studies JMJHlM HE HE HJ HI HE 

in hi he m he he nr ErEr, to which we shan 

refer in the following. Those studies are generally two- 
dimensional (except m HE HE]), with a prescribed ve- 
locity (except |44| where the shear stress is prescribed), 
with a fixed volume (except |H HE HI HE Eg where the 
pressure is prescribed), and between rough walls (except 
HH HE HE HE 113 where periodic boundary conditions 
are applied along the transverse direction). 



B. Contact law 

Let us consider two grains i and j in contact, of di- 
ameter d^.j, mass rriij, centered at position f^j, with 
velocity v^j and rotation rate ui^j. Let Hij denote the 
normal unit vector, pointing from i to j (Hij — rij/\\vij\\ 
with the notation Hij = Hj — Hi), and ty a unit tan- 
gential vector such that (Hij,tij) is positively oriented. 
We call Fij — NijHij + Tijtij the contact force exerted 
on the grain j by the grain i. The contact law re- 
lates the normal, Nij, and tangential, T^, components 
of the contact force to the corresponding components 
of relative displacements and/or velocities. The rela- 
tive velocity at the contact point is equal to Vij = 
Hi — Hj + l/2(diUJi + djUjj)tij. Its normal component 



N?j = k n hij. (4) 

Eqn. (0} involves a constant normal stiffness coefficient 
k n , the value of which is independent of disk radii. 

The normal viscous force opposes the relative ap- 
proaching or receding velocity : 

Nij = Qjhij. (5) 

Qj is related to the normal restitution coefficient e in a 
binary collision, and chosen such that e is constant for 
all contacting pairs, whence 

-21ne jrriij 

Qj= / a ■ , 2 =V IT ■ (6) 
V 7T 2 + hr e V fc » 

where = rriimj / '(m, + rrij). 

Physically, Eqn. (@J can be regarded as a simplified 
version of the Hertz law > N e oc h 3 / 2 , while the viscous 
dissipation mi ght stem from the visco-elasticity of the 
grain material [44( • 

The total normal force might be either repulsive or at- 
tractive, due to the viscous contribution. We could check 
that setting Nij to zero whenever it becomes attractive 
(Nij < 0) has but a negligible effect on the simulation 
results l45l. 



2. Friction 

The Coulomb condition in the contacts involves the 
coefficient of friction between grains /i, and is enforced 
with the sole elastic part of the normal force : 

\Tij\Kfim, (7) 
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polydispersity fi e kt/k n k 
~ ±20% 0-0.8 0.1-0.9 0.5 1(F 

TABLE II: List of material parameters. 

To this end, the tangential component of the contact 
force is related to the elastic part 5fj of the relative tan- 
gential displacement Sij , 

Ta = hStj, (8) 

with a tangential stiffness coefficient k t . 8fj is defined 
by : 

d 3_ = f if \T %3 1 = M iV?. and T^V? > 
dt VVg otherwise 

The contact is termed "sliding" in the first case in 
Eqn. © (the condition that Ty and V? share the same 
sign ensuring a positive dissipation due to friction) and 
"rolling" in the second case, kt is of the same order of 
magnitude as L As it has a very small influence 

on the results [39], it was fixed to k n /2 in all our calcu- 
lations. 

Table [H] gives the list of material parameters (antici- 
pating on Sec. Ifffl we indicate the value of the dimen- 
sionless parameter k). Most simulations were done with 
e = 0.1 and /i = 0.4, which implies a rather strongly 
dissipative material, favors dense flows, and seems fairly 
realistic. Other friction coefficients and restitution coef- 
ficients were also studied. 



C. Simulation method 

The interaction law being chosen, numerical simu- 
lations are carried out with the molecular dynamics 
method, as in refs. l47l |48[ . The equations of mo- 
tion are discretized using a standard procedure (Gear's 
order three predictor-corrector algorithm ^j). The time 
step is a small fraction (1/100) of the duration r c of 
a binary collision between two grains of mass m (r c = 

v /m( 7 r 2 + ln 2 e)/(4fc„)). 



of the simulated system. If the length and mass scales 
d and m are obvious, there are three candidates for the 
time scale : the shear time I/7, the inertial time y/m/P 
(that is to say the characteristic displacement time of a 
grain of mass m submitted to a pressure P), and the 
collision time r c . However, dimensional analysis predicts 
that the behavior will only depend on six dimensionless 
numbers. Apart from /z, e and k t /k n , we propose the fol- 
lowing choice of the three other dimensionless numbers: 

i) The first g p /\/mk n is associated to the normal mo- 
tion of the wall controlling the pressure. A small 
value signifies that the time scale of the fluctua- 
tions of H is imposed by the material rather than 
the wall, and that the wall "glues" to the material. 

ii) The second k = k n /P is associated to the rigidity 
of the grains. It is inversely proportional to the 
normal deflection h of the contacts for a confining 
pressure P. 

iii) The third (V/H)y/m/P describes the shear state, 
through a combination of the three global parame- 
ters (velocity V, pressure P and height H), which 
means that it is not necessary to vary indepen- 
dently those three parameters. 

The first dimensionless number is equal to 1 in our sim- 
ulations. We have not studied its influence, but consider 
that is is small in this range. We shall not refer to it in 
the following. 

Most of our calculations are restricted here to the limit 
of rigid grains 01 j K ~~ * 00 j as we chose k = 10 4 (unless 
specified otherwise, as in Sec. HVD~2|) . From studies in 
the quasi-static regime, this value is known to be large 
enough for the coordination number to show little varia- 
tion as k is further increased (see, however, Sec. lVII A")) . 
With k = 10 4 , the typical ratio h/d oc k -1 is so small 
that elastic deflections h stay negligible in comparison 
with the gaps between neighboring grain surfaces that 
determine the amplitude of rearrangement events 0] (a 
more stringent condition than h <C d). In practice, a 
suitable definition of k for particles with Hertzian con- 
tacts (such that h/d (x kT 1 with a coefficient of order 1) 
is k = (E/P) 2 / 3 , where E is the Young modulus of the 
material the grains are made of. For glass (E — 70 GPa) 
and P — 10k Pa (the pressure due to the weight of a 50 
cm thick layer) one has k ~ 37000. 



III. DIMENSIONAL ANALYSIS 

In order to analyze the results, we recall the list of 
parameters describing the material and the shear state. 
The grains are described by their size d, mass m, stiff- 
ness parameters k n and k t , and coefficients of restitution 
e and of friction fj,. The shear state is described by the 
prescribed pressure P and the average shear rate V/H, 
and by the viscous damping parameter g p . In the fol- 
lowing, we shall express the results in the natural units 



A. Inertial number I 

In the following, we study in detail the influence of 
the third dimensionless number. It is a global number 
I g , at the scale H of the sheared layer. Replacing the 
global shear rate V/H by the local one 7 defines the local 
analogous quantity : 
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As the ratio of inertial to shear times, I measures the 
inertial effects, and will be called inertial number in the 
following. This number already appeared in the colli- 
sional regime (Sec. II Bf) . In a three dimensional situa- 
tion, the inertial number is equal to 7y m/Pd. We no- 
tice that, introducing the mass density p g of the grains, 
the definition / = jdw p g /P does not depend on the di- 
mensionality of the system. 

In a homogeneous system, without sliding velocity at 
the wall, both definitions, global and local, are equiva- 
lent. But in a heterogeneous system, it is necessary to 
take into account the variations of the local inertial num- 
ber along y: I(y) = j(y)y/m/P(y). 

We are going to show that this dimensionless number 
might well be the fundamental quantity to describe the 
rheology of granular materials. First of all, it is the cor- 
rect quantity to characterize the flow regime. Thus, the 
terminology of "slow" versus "rapid" granular flows is not 
correct, since it is necessary to combine shear rate and 
pressure to characterize the inertial effects. A small value 
of I (small 7 and/or large P) corresponds to a regime 
where the grain inertia is not relevant : this is the "quasi- 
static" regime (Fig. H( a ))- Inversely, a large value of / 
(small P and/or large 7) corresponds to the "inertial" 
or "dynamical" regime, which may be described by the 
kinetic theory (Fig.[T](b)). Varying / allows to study the 
progressive transition between those two regimes. The 
range of inertial number which we have studied goes from 
6 • 1CT 4 to 0.3. 



B. Comments 

It might be pointed out that the inertial number is di- 
rectly related to the Bagnold number Ba, which can be 
defined as the ratio of the typical kinetic energy m(dj) 2 
of one grain to the characteristic frictional dissipation 
fiPd 2 0. Some authors [28ll39| chose to control the solid 
fraction v rather than the pressure, and therefore, rather 
than I and k, used the pair of dimensionless numbers v 
and a — 7/ yjk n /m = I/^/k, as variables characteriz- 
ing the state of the granular material in steady homoge- 
neous shear flow. The latter dimensionless combination 
a may be viewed as the ratio of the collision time to 
the shearing time, or as the shearing velocity divided by 
the sound velocity (Mach number 3!J|). Both choices are 
perfectly legitimate, as dimensional analysis predicts ei- 
ther S/P = /i(7, k) and v = f 2 (I, k), or S/P = f 3 (u, a) 
and P/kn — fi{v,a), both results being equally valid. 
The choice of / and k can however be deemed more con- 
venient for several reasons. First, the variation of the 
results with v, regarded as a control parameter, is ex- 
tremely fast. Each material possesses a critical packing 
fraction v c , in the sense of Sec. II Al above which it does 
not flow, unless stresses are so large that the elastic com- 
pression of contacts compensates for the difference v — v c . 
Below v c , on the other hand, a continuously sheared gran- 
ular system is free to flow with a negligible shear stress, 



unless the velocity is high enough to build a significant 
pressure. One should therefore monitor v with great ac- 
curacy to observe ordinary stress levels. This renders the 
comparisons between different granular systems difficult, 
as one would need to know in advance the value of the 
critical density for each of them. Furthermore, the limit 
of rigid grains becomes singular, as all values of v above 
v c are strictly forbidden for 7 7^ 0, while the properties 
of shear flows with v < v c simply scale with 7 in that 
limit 01 ■ Conversely, if one uses I and k as control pa- 
rameters, no singularity enters any of the relevant results 
in the I — » or k — > 00 limits, and different materials 
should exhibit similar (if not quantitatively identical) be- 
haviors for the same values of these parameters (which 
thus define roughly "corresponding states"). It should 
also be pointed out that experimental conditions usually 
determine stress levels, rather than densities. 



IV. HOMOGENEOUS SHEAR STATE 

A. Preparation of steady shear states 

The first kind of preparation (which has been used 
most of the time) consists in starting from an initial con- 
figuration where the disks are randomly deposited with- 
out contact and at rest between the two distant walls, 
which provides an average solid fraction of 0.5, and then 
in applying the pressure to the wall while slowly shearing 
the granular material. When the pressure on the walls 
reaches the prescribed value, the prescribed velocity is 
applied. 

The second kind of preparation consists in starting 
from a very high solid fraction (of the order of 0.8), ob- 
tained by a random deposit followed by a cyclic com- 
paction with frictionless grains, and then to introduce 
friction between grains and to start the shear with the 
prescribed velocity and pressure. 

The two kinds of preparation allow to start either from 
a loose state (the first case), or from a dense state (the 
second case). Occasionally, we have used a third kind 
of preparation consisting in starting from localized shear 
states near one of the walls (obtained by applying gravity, 
see Sec. lIXf) . 

We then look for a steady flow, characterized by con- 
stant time-averaged quantities of the flowing layer, like 
kinetic energy and solid fraction. We have observed that, 
after a sufficient amount of time, the three kinds of prepa- 
ration lead to the same shear state. We deduce that there 
is no influence of the preparation on the steady flow char- 
acteristics. 

All the simulations converge to an average steady state. 
But the relative fluctuations of the measured quantities 
can be very different. In the following, we consider that 
the shear state is continuous if they are smaller than 10%. 
Otherwise, the flow is called intermittent. This happens 
in the quasi-static regime, for I < 0.001 (see Sec. I VI All . 
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When a steady state is reached, the simulation is car- 
ried on during a sufficient amount of time At, so that 
the typical relative displacement of two neighboring lay- 
ers is larger than ten grains (7 At > 10). In this steady 
state, we consider that the statistical distribution of the 
quantities of interest (structure, velocities, forces. . . ) are 
independent of t and x, so that we average both in space 
(along x) an in time (considering 200 time steps dis- 
tributed over the period At). 

B. Profiles 

We now show that the granular material is completely 
sheared and that the shear is homogeneous. To this 
end, we analyze the profiles of solid fraction, velocity and 
stress components. Typical profiles are shown in Fig. [21 

1. Solid fraction and velocity 

The profiles of solid fraction v(y) and shear rate 7(2/) 
are strongly oscillating near the walls, but we may con- 
sider them as approximately constant in the central part 
of the sheared layer (Fig. El (a) and (b)). This allows to 
define average solid fraction v and shear rate 7. 

We notice that those profiles obtained for controlled 
pressure are in agreement with the one measured at fixed 
volume [50| . 

We have also measured the profile of average rotation 
rate ui(y) (Fig. (b)), and observed, as in other flows 
and quasi-static deformations EH EE EE E~J E3- 
the following relation ui(y) — —57(2/). 

2. Stress tensor 



The stress tensor E is the sum of two contributions 




The first term ("contact"), usual in static of gran- 
ular materials, is associated to contact forces between 
grains [5?l l5Sj . The second term ("fluctuations"), usual 
in fluid mechanics (Reynolds tensor), is associated to the 
velocity fluctuations of the grains Svi : 

(12) 

£ = XT? Th=x m i^i ® 

A third contribution, associated to the rotation of the 
grains, has been introduced in 15911 . We have observed 
that is is usually insignificant (2y|, so that we shall not 
discuss it in the following. We discuss here the total 
contribution and shall analyze the two contributions in 
Sec. El 




y/(v/H), -2co/(v/h) 




S(y)/P, P(y)/P 



FIG. 2: (Color online) Homogeneous shear state : (a) Solid 
fraction v(y), (b) Shear rate 7(3/) (—) and rotation velocity 
—2u{y) (o), (c) Pressure P(y) (—) and shear stress S(y) (o) 
(V = 1,P = 25, H/d w 40, e = 0.1, (j, = 0.4;. 



Using the same averaging methods than for solid frac- 



E = E C + S / . (11) y/d20 
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tion and velocities 55], it is possible to calculate the 
profiles of the different components of the stress ten- 
sor. We observe that E™, ~ E yy , like for dense flows 
down inclined planes [42, H3, [55| . Consequently the pres- 
sure P = (Yi X x + Ej /a )/2 ~ ^yy In the following we 
call S = -T, xy . The Fig. (c) shows the profiles P(y) 
and S(y). Those components arc approximately constant 
along y, as expected. 



3. Inertia! number 

From the profiles of shear rate and pressure inside the 
flow, we deduce the profile of inertial number I(y). ft 
is approximately constant in the center of the sheared 
layer, so that it is possible to define an average 7 (in 
Fig. |1 7 w 0.004). 



C. Behavior near the walls 

Particular behaviors of solid fraction, velocity and 
stress components are observed near the walls. First, 
the walls induce a structuration of the granular mate- 
rial in the 5 first layers, evidenced by the oscillations of 
the solid fraction and shear rate. Second, the frustra- 
tion of the rotation of the flowing grains in contact with 
the rotationless glued grains of the roughness induces a 
sudden variation of the average rotation. This may be re- 
sponsible for the difference of normal stress components 
{Y, xx ^ Eyy) in the first two layers (not shown). Third 
the granular material is agitated by the roughness, which 
increases the shear rate (see Sec. IVIJI . Consequently, all 
the average quantities are measured in the central part 
of the sheared layer, excluding the 5 first layers near the 
walls and 77 is chosen large enough so as to limit those 
wall effects. 

However, we have tested the homogeneity of the shear 
in the case of a thin layer (77/<i sa 5) [26J. Then the 
granular material is structured on its whole width and 
completely sheared. We shall come back to this case in 
the following. 



1, Influence of the polydispersity 

In the case of a very small polydispersity (< 1%), we 
have observed [2(| that the granular material crystallizes 
near the walls. The shear zone then reduces to a ten 
diameters thick central layer, where the velocity profile 
is linear. When the polydispersity is increased, the shear 
zone extends first to the moving wall (5%) then to the 
whole layer (10%). 



2. Influence of the rigidity 

We expect that the correlation length of the strain field 
becomes smaller when the grains become softer. This 
might be responsible for a localization of the shear near 
the moving wall. In order to test this idea, we have stud- 
ied the influence of the dimensionless rigidity number k 
in a large system (H/d = 100) for I g = 2 • 10~ 4 . For 
k = 1 000, we observe homogeneous shear states, whereas 
for k — 40, the shear becomes localized |26[ . 



V. CONSTITUTIVE LAW 

Those homogeneous states are a great advantage to 
measure the constitutive law of dense granular flows. 
This would be much more difficult in other flow geome- 
tries (annular shear, vertical chute, inclined plane, heap- 
flow. . . ) [2{| where the flow is heterogeneous, so that the 
constitutive law cannot be identified on the whole flow- 
ing layer, but should be studied locally, in the presence 
of gradients. 

In those homogeneous states where the inertial num- 
ber 7 is prescribed, the solid fraction and the shear stress 
adjust in response to 7. Consequently, we shall first mea- 
sure the dependences of two fundamental dimensionless 
quantities, the solid fraction and the effective friction co- 
efficient, as functions of 7. Then we will study the in- 
fluence of the mechanical properties of the grains, the 
coefficients of restitution e and friction ji. 



D. Limits of homogeneity 

Even when starting from a localized velocity profile, we 
have observed a relaxation toward a homogeneous shear 
state. This observation is in contrast with other studies, 
where a shear localization is observed [53, [35|, [3?1 l38| . 
In our simulations, we have observed signs of localiza- 
tion in the quasi-static regime where the flow becomes 
intermittent fSec. IVl|) . and once the stress distribution 
becomes heterogeneous (SeclIXjl. We now briefly discuss 
two other factors of localization, which are associated to 
the granular material itself, the monodispersity and the 
softness of the grains. 



A. Dilatancy law 

We call "dilatancy law" the variations of the average 
solid fraction v as a function of the inertial number 7 
(Fig. EJ. We observe that v decreases approximately lin- 
early with 7, starting from a maximum value v max : 

~ v max - al, (13) 

with v max ~ 0.81 and a ~ 0.3 (for // = 0.4). The er- 
ror bar (independent of 7) corresponds to the statistical 
dispersion inside the layer. 
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FIG. 3: Dilatancy law : (a) Linear scale (b) Logarithmic 
scale (/j, = 0.4, various e and k). 



cussion of the influence of the roughness on the flow (in- 
clined plane and vertical chute), we refer to l55l l60 | . 
In the following we shall only discuss the effective friction 
coefficient in the volume of the flowing layer. 

We call "friction law" the variations of the effective 
friction coefficient fi* (averaged over the width in the cen- 
tral part of the flowing layer) as a function of I (Fig. 0}. 
We observe that [i* increases approximately linearly with 
J, starting from a minimum value ix* min ■ 



/'r, 



bL 



(14) 



with \i* min — 0.22 and b ~ 1.0 (for /i == 0). The error bars 
(independent of /) correspond to the statistical disper- 
sion inside the layer. We also observe that \i* tends to 
saturate for / > 0.2. Within the error bars, it is difficult 
to be more precise about those dependencies. A more 
careful measurement is deferred for future work. 
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B. Friction law 

The effective friction coefficient has been defined as 
the ratio of the shear stress to the pressure inside the 
material //* = S/P. 

It could also be defined as the ratio of the (total) tan- 
gential and normal forces on the wall /it* = T/N. We 
have observed [2fJ that aC is slightly larger than /i*. 
Some simulations have been carried out on to test the 
influence of the roughness, by taking glued grain on the 
wall twice as small (R = 0.5) or twice as large (R = 2) as 
the flowing grains, for the same I g . This size ratio has an 
influence on the sliding velocity at the wall : it becomes 
noticeable for R = 0.5 and decreases when R increases, 
since the grains close to the walls are trapped by the 
roughness. However, at distance from the walls, the flow 
remains homogeneous, but the shear rate, and hence J, 
decreases when R decreases. Furthermore, the effective 
friction at the wall decreases when R decreases. All in all, 
aCCO seems independent of R. For a more detailed dis- 




FIG. 4: (Color online) Friction law : (a) Linear scale (b) 
Logarithmic scale (fj, ^ 0, various e and k). 



We now compare this friction law with other works. 
We first notice that the increase of n* with / is contrary 
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to the well documented decrease of the friction coeffi- 
cient with the velocity in the quasi-static regime |6lj . 
However, in those studies, this softening is interpreted 
as a consequence of the renewal of the population of 
asperities at a microscopic scale, or as an effect of hu- 
midity 62] . Those effects do not come into play in our 
study. As a matter of fact, this friction law was already 
observed in previous discrete simulations |3l|. and par- 
tial observations (experimental or numerical) of the vari- 
ation of /i* with the shear rate, the pressure or the solid 
fraction, consistent with our observations, may be found 

in^lHimmiHililil In terestingly, the 
inclined plane geometry allows to prescribe both the ef- 
fective friction and the pressure, through the inclination 
9 of the plane and the height H of the flowing layer. 
Consequently, the measure of the superficial velocity V 
as a function of these two parameters provides a measure 
of the friction coefficient at the base as a function of I g 
(which is proportional to V/H 3 / 2 ) |69ll70| . Those obser- 
vations are in good agreement with the previous friction 

law mEt 



C. Comments 

As a conclusion, in the limit of rigid grains, the classi- 
fication of the flow regimes depends on the single dimen- 
sionless number I. In the quasi-static regime (/ < 10~ 2 ), 
the granular material is very dense, close to the maximum 
solid fraction v max , and the effective friction coefficient is 
close to its minimum value [i* min - In the dynamic regime 
(/ > 0.2), the dilatancy becomes strong and the effec- 
tive friction coefficient seems to saturate. The transition 
between those two regimes is progressive. In the interme- 
diate regime (10~ 2 < / < 0.2), we observe approximately 
linear variations of the solid fraction and of the effective 
friction coefficient as a function of / (Eqs. and 114fl ). 

In the case of a thin layer (see Sec. llV Ofr . the dilatancy 
and friction laws are not affected. The parameters remain 
the same, except for the solid fraction which is smaller 
than for a thick layer : v max — 0.82 instead of 0.84 for 
[i = and v m ax — 0.78 instead of 0.81 for fi — 0.4. 

In the case of a very small polydispersity (see 
Sec. II V Pfl . the friction law is preserved in the central 
sheared layer. For larger polydispersity (< 50%), the 
dilatancy and friction laws are not affected. 

When taking into account the elasticity of the grains, it 
is natural to draw a diagram of the flow regimes as a func- 
tion of the two dimensionless numbers, already quoted in 
Sec. Ens a = I/s/k and v [HHl]. This leads to iden- 
tify three regimes : elastic quasi-static, purely inertial 
and elastic-inertial. This last regime corresponds to very 
soft grains (k < 100) and is not accessible in our study, 
where we stay in the limit of rigid grains. 

Furthermore, various studies where the volume rather 
than the pressure was prescribed [2^, 0, |3j| |§3 have 
evidenced a transition between the quasi-static and the 
inertial regimes around a critical solid fraction. In our 



study where the pressure is prescribed, the solid fraction 
adjusts to the inertial number /, so that the transition is 
not accessible. 

D. Constitutive law 

We shall now show how these dilatancy and friction 
laws allow to deduce the constitutive law of the material, 
that is to say the dependencies of the pressure and shear 
stress on the shear rate and solid fraction : P{v, 7) and 
5(i/, 7). 

1. Pressure 

From the definition of I (Eqn. (JTJJJ) and the dilatancy 
law (Eqn. (|13|) ). the pressure may be expressed as a func- 
tion of the shear rate and the solid fraction : 




V - V 

max 



FIG. 5: -P/7 2 as a function of (v max — v). The straight line 
indicates a slope of —2 (fj, — 0.4, various e). 



2. Shear stress 

From the definition of I (Eqn. I|10(l ) and the friction 
law (Eqn. JHJ), the shear stress may be expressed as a 
function of the shear rate and the pressure : 

S(P,i) = fx* min P + bV^Pj. (16) 

Using the previous expression of P (Eqn. I|15[)). it is 
also possible to express the shear stress as a function of 
the shear rate and the solid fraction : 
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S(P,u,i) = ^ mln P+- ( -rmf, (17) 

or, eliminating P : 

^7)=^£^F^ 2 > ( 18 ) 

with the solid fraction v* — v max + a^* min /b [y* w 0.86 
for fj, = 0.4). 



3. Comments 

We first notice that the dilatancy and friction laws, 
measured in the whole range of regimes, from the quasi- 
static to the dynamic, make the link between the known 
results in the two extreme regimes recalled in Sec. |I] In 
the quasi-static regime, [i* min and v max may be identi- 
fied with the internal friction tan <f> and the critical solid 
fraction v c in the critical state (see Sec II A\i . 

In the intermediate regime (I < 0.2), the expression of 
the stress components is analogous to the expression for a 
homogeneous system in the collisional regime (Eqs. ©)■ 
The dependence on the square of the shear rate, similar 
to the original conclusions of Bagnold for concentrated 
sus pen sions |72j . is a consequence of dimensional analy- 
sis • The dependences on solid fraction are described 
by the following dimensionless functions : 




(19) 



We recall (see Sec. II B|l that the kinetic theory predicts 
a divergency in 1/(1 — v) 2 from the asymptotic behavior 
of the pair correlation function. However, when v > 0.67 
(so-called gel transition), then starts a regime of multiple 
collisions, strongly correlated, and the divergency rather 
seems in l/(v max — v) |4jj. When the material is 
sheared, under the effect of cooperative rearrangements 
of caged grains, an even stronger divergency of the viscos- 
ity has been conjectured It seems that the precise 
form of these divergencies are decisive to describe the 
shape of the velocity profiles and for the jamming pro- 
cess jT^, 0| . Our quantitative determination is then a 
precious information for the modelling of the dense gran- 
ular flows. 

We also notice that the expression l|16|) of S corre- 
sponds to a visco-plastic constitutive law, similar to the 
"frictional-collisional" decomposition of the stress tensor, 
with a contribution associated to maintained contacts, 
and a contribution associated to collisions |4jJ, [fH 
FtI l7S| . In the viscous term, we notice that the appar- 
ent viscosity b\JmP is proportional to the square root of 



the pressure. The interpretation it that the typical mo- 
mentum rwyd is exchanged with the inertial time scale 
yJm/P over a surface of the order d 2 . 

We think that the formulation of the constitutive law 
through the dilatancy and friction laws is simpler to use, 
since it avoids the treatment of divergency near jamming, 
which might be a problem in fluid mechanical numerical 
simulations. 



E. Influence of the mechanical parameters 

We shall now describe the influence of the mechanical 
properties of the grains (fc„, \x and e) on the dilatancy 
and friction laws. 



1, Influence of the elasticity of the grains 

In the limit of rigid grains, that is to say for k > 10 4 , 
we have not observed any influence of the stiffness coef- 
ficient k n on the constitutive law. There is however an 
influence on the coordination number (see Sec. IVII Afl . 
Furthermore, we have shown that softer grains favor lo- 
calization of the shear (see Sec. HVD~2"|) . 

2. Influence of the local friction coefficient 

The local friction coefficient // has a significant influ- 
ence on the dilatancy law. The solid fraction remains 
a linearly decreasing function of I, but both parameters 
v max and a depend on [i. 

The variation of a is not simple: we measure a = 0.38 
for n = 0, 0.31 for fi = 0.4 and 0.37 for /Lt = 0.8 H|. Fig.© 
(a) indicates that v max is a decreasing function of fi (it is 
not purely geometrical in the quasi-static regime). This 
shows that the solid fraction, from the critical state to the 
collisional regime, depends on the frictional properties of 
the material. We notice that our measurements are in 
agreement with other observations 0, |3^, . 

The influence of fi on the friction law is less significant, 
except for frictionless grains (fi = 0) . The Fig. where 
fi varies between 0.1 and 0.8, shows that fi has nearly no 
influence on ji* in this range. This variation is more sig- 
nificant for small /. The Fig. (b) shows more precisely 
the variation of the effective friction as a function of /i 
in the quasi-static regime. There is strong variation be- 
tween fj, = and jU = 0.4, but above \i = 0.4 the effective 
friction remains constant. 

As a conclusion, in the case of frictionless grains, the 
friction law keeps the same tendency but is shifted to- 
ward smaller values of friction (Fig. |7J). However, we 
observe a saturation for / > 0.1 if e = 0.9, which will be 
discussed in Sec. IV E4l For / < 0.1, the linear approxi- 
mation (Eqn. (|14|) 1 is in trouble : it is rather a sub- linear 
dependency, and \x* min — 0.11. So the case of frictionless 
grains is singular and deserves a specific treatment. 
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FIG. 6: Influence of fj, on the critical state (I = 4 • 10~ 3 and 
e = 0-1): (a) Maximum packing fraction v max , (b) Effective 
friction coefficient [i* m i n , (c) Coordination number Z max . 



Starting from both variations of solid fraction and ef- 
fective friction as a function of the local friction coef- 
ficient fx, it is tempting to draw the variations of the 
effective friction as a function of solid fraction instead of 
inertial number. This is done on Fig. |SJ As a matter of 



FIG. 7: Influence of n on the friction law (\i — (U), fi ^ 
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FIG. 8: (Color online) Variation of the effective friction co- 
efficient as a function of solid fraction. 



fact, Eqs. (fT5)l and fl%|) predict 



H*{y) = -{v* - is), 
a 



(20) 



where v* was previously defined. This new represen- 
tation of the results evidences a collapse of the data 
(even if a, b, n* min and v max vary separately with /z, 
v* = Vmax + seems approximately constant). It 

appears that /i* becomes nearly independent of fj,. This 
master curve is made of complementary zones of high 
solid fraction for frictionless grains, and smaller solid 
fraction for frictional grains. It is noteworthy that a 
small variation of solid fraction (of the order of 10%) 
is enough to induce a variation of effective friction by a 
factor 4 ! This decrease of the effective friction when 
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the solid fraction increases is contrary to the observation 
in static compaction [8(i|. but was previously observed 
under shear |6fij . 

3. Influence of the restitution coefficient 

The Fig.0]shows that there is no influence of the resti- 
tution coefficient e for frictional grains. For frictionlcss 
grains, the comparison between slightly (e = 0.1) and 
strongly (e = 0.1) dissipative grains reveals (Fig.[5J that 
there is an influence, limited to the dynamic regime, for 
I > 0.1. Then, as the dissipation decreases, the dilatancy 
is less pronounced and the effective friction saturates (see 
also |33p. 



elastic collisions (then the average contact time tends to 
the collision time |H3)- I n ^ ms dense limit, the Eqs. © 
predict a value of the effective friction independent of I 
but dependent on e, from the values of the pre-factors 
Ai [3 : 

*® = —zr = i\r*r) Vl - e - (21) 

For e = 0.9, this predicts fi* — 0.29 which is in fairly 
good agreement with the value measured for / = 0.2 
(0.26). The small difference between those two values 
may be due to the influence of the walls which induce a 
sliding velocity and gradients of the fluctuations |8l|. 
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VI. FLUCTUATIONS 

Up to now, all the quantities we have studied (solid 
fraction, velocities, forces) have been averaged in space 
and time. In fact, they are heterogeneous in space and 
fluctuate in time. The fluctuations of the forces have 
already been thoroughly studied [H HHH H2| . We 
have observed [2(| that the fluctuations of the forces are 
of the order of a few percents in the quasi-static regime, 
increase up to a factor 3 to 4 in the dynamic regime, and 
that they are larger near the walls than inside the sheared 
layer. We shall now discuss in detail the fluctuations of 
the motion of the grains. 




0.20 



FIG. 9: Influence of e for /i = (e = 0.9 (M), e = 0.1 (o)) 
(a) Dilatancy law, (b) Friction law. 



4- Collisional limit 

The situation (I > 0.1, /i = and e = 0.9) corresponds 
to the dense limit of the kinetic theory, with binary quasi- 



A. Intermittencies in the quasi-static regime 

Various studies have shown that granular flows become 
unstable in the quasi-static regime. When the velocity 
is prescribed, one goes from a continuous flow regime to 
stick-slip [62ll83ll84 |. When the shear stress is prescribed, 
one observes an hysteretic and abrupt fluid-solid transi- 
tion [8^ . We also notice that the observation of shear lo- 
calization in discrete simulations of plane shear without 
gravity corresponds to the quasi-static regime [35L l38| . 
Space-time correlations of the motion of the grains have 
been observed |H HI |H HH HU . Several rheo- 
logical models have been prop osed to describe those ob- 
servations m m m m hii m m m @ m : trans- 
mission of forces at the scale of correlated clusters, two- 
phase fluid model with order parameter, activation of 
rearrangements through the fluctuations of velocity or 
forces. 

Our study (see j2(J for more details) shows that the 
flow is continuous in the intermediate regime, but be- 
comes intermittent in the quasi-static regime (for / < 
I ~ 0.003). Then the time averaged shear is homoge- 
neous, but the instantaneous velocity profiles show that 
the layer oscillates between two states : a first one where 
the whole layer is at rest, caught to the fixed wall, except 
for a thin layer swept along by the moving wall, and a 
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second one where the whole layer is swept along by the 
moving wall except for a thin layer caught to the fixed 
wall. Those two extreme states have a very short dura- 
tion, and the system is most of the time in an intermedi- 
ate state, where the shear is approximately homogeneous 
in the whole layer. The pictures of the fields of velocity 
and velocity fluctuations show that the spatial correla- 
tions have a vortex like shape with a size between 3 and 
4 grain diameter and are of a short duration. Conse- 
quently, the total kinetic energy fluctuates in time, with 
sudden peaks associated to the two extremal states, and 
slow variations associated to the intermediate situations. 
The Fig. ^] indicate that the relative fluctuations of the 
kinetic energy increase when / decreases and saturate in 
the quasi-static regime, according to : 



<E C > 
with I ~ 0.003. 



£ forK/o, 
% for I > In, 



(22) 
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the motions of the grains are strongly correlated [42j. 
The definition itself of those fluctuations is a problem, 
since it has been shown that they depend on the av- 
eraging time scale in the quasi-static regime [2j|. Our 
analysis (long time scale) takes into account both the 
small fluctuations around the mean motion (in the "cage" 
formed by the nearest neighbors 96] ), and the large fluc- 
tuations associated to collective motions in the quasi- 
static regime. 

We notice that the quantities defined by Eqs. ill'MII are 
local, i.e. function of y. Those velocity fluctuations are 
uniform at the center of the sheared layer, but increase 
near the wall [2EIH3- According to various models, this 
agitation of the granular material near the rough walls 
is responsible for the increase of the shear rate (observed 
in Fig. H (b)) HHH El HI. In the following, we 
consider the average velocity fluctuations in the central 
part of the sheared layer, excluding the 5 first layers near 
the walls. 
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FIG. 10: Relative fluctuations of the kinetic energy as a func- 
tion of I . Influence of H . 



B. Velocity fluctuations 

We now discuss the fluctuations of the translation ve- 
locity Sv and of the rotation velocity 8u> of the grains, 
which measure the agitation of the granular material. 
They are defined by : 



§©/y 




6v{y) = sj < v(y) 2 > - < v(y) > 2 , 
Su(y) = V< ^(y) 2 > - < uj(y) > 2 . 



(23) 



They are a natural variable in the description of col- 
lisional granular flows (see Sec. II B(l . Their influence is 
more difficult to analyze in the case of dense flows, where 



FIG. 11: Relative velocity fluctuations as a function of I : 
(a) Translation velocity, (b) Rotation velocity (various param- 
eters). 



Dimensional analysis suggests to analyze the variations 
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of the dimensionless quantities Sv / jd and Sui/'j as func- 
tions of /. We consider -yd as the natural scale of transla- 
tion velocity, and 7 as the natural scale of rotation veloc- 
ity (since w = — ^7, see Sec. II V Bill . Those variations, 
drawn on Fig. 1111 evidence two scaling laws, independent 
of the parameters of the system : 



8v ^ 1 T—ct 
id — 3 ' 



5u ^ j-j3 

7 ' 



(24) 



with a ~ \ and (3 ~ -|, to which we shall refer as trans- 
lation and rotation scaling laws. 

The velocity fluctuations are significant (larger than 
1) in the intermediate regime and very significant (larger 
than 10) in the quasi-static regime. In comparison, ki- 
netic theory (see Eqs. ||TJ and predicts that Sv/jd — 

yjA s /Ar that is to say 32 (i^e^) — 1<3 > for e = °- 9 ' 
which corresponds to the order of magnitude which is 
measured in the dynamic regime. 

Large values of Suj/j have been observed experimen- 
tally in the quasi-static regime [HJ and in dense flows 
down inclined planes l55ll. and might be due to the frus- 
tration of the rotation |98| . 

We give two interpretations of the value of the expo- 
nent a. The first explanation [2^] consists in analyzing 
the motion of one grain as a succession of shear phases 
of duration 1/7 with a velocity jd and of sudden rear- 
rangements with a velocity d\[Pjra of duration ^JmjP . 
This leads to 4^ ~ I^ X ^ 2 jlfj- The second explanation 
relies on an energetic argument. In homogeneous shear, 
the work of the shear stress is balanced by the dissipation 
rate S7 = Y. If V describes the dissipation of the fluc- 
tuating kinetic energy mSv 2 /2 during the inertial time 
yJm/P, we obtain : Sv/jd ~ yj2p*(I)I~ 1/2 . For both 
interpretations, the order of magnitude of the pre-factor 
is consistent with the observation. 

We now show that the translation scaling law is con- 
sistent with the variations of the relative fluctuations of 
the kinetic energy (Eqn. I|22|l ). Let us call p = p g v the 
average solid fraction of the granular materials and (E c ) 
the average kinetic energy by unit length, which is domi- 
nated by the translational part : (E c ) = f J^ 1 ( / jy) 2 dy ~ 
(pV 2 H). In the quasi-static regime, where the system os- 
cillates between two localized flows, AE C ~ (pV 2 H)/2, 
so that AE C /(E C ) ~ 1. In the dynamic regime, AE C ~ 
2 Jo 3 v2 dy- Using the translation scaling law for this 
last quantity, we get : AE C ~ pV 2 /(HI), so that 
AE C /(E C ) ~ 1/(H 2 I). This is in agreement with the 
dependencies on / and H observed in Fig. UHl 

Furthermore, the translation scaling law provides an 
estimation of the Reynolds contribution to the stress ten- 
sor (see Eqn. (0): T, f /T, ~ (2v/9w)I. This shows that 
in the intermediate regime (/ < 0.1), this contribution 
remains smaller than 1%, so that the contribution of the 
contact forces S c remains dominant. 



C. Consequences for the constitutive law 

The translation scaling law may also be written : 
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7 X / 2 (p/ m )i/4. 



(25) 



Consequently, when the pressure is prescribed, the ve- 
locity fluctuations vary like 7 1 / 2 , instead of 7. In the 
annular shear geometry, the pressure is constant along 
the radial direction, and an exponent close to 1/2 has 
been measured experimentally 18(| . 

If we introduce the velocity fluctuations in the consti- 
tutive law, like in the collisional regime (see Eqn. Q), 
we obtain : 
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(26) 



Within this formulation, we notice a stronger diver- 
gency of the viscosity near the maximum solid fraction, 
like in the model inspired by the glassy dynamics [74|. 



VII. CONTACT NETWORK 

In addition to the macroscopic quantities which have 
been discussed, the discrete simulations provide informa- 
tion on the contact network. Its strongly heterogeneous 
character both in s pac e and time has already been dis- 
cussed in detail pi fsl Isl Isl lil Isl Isl Isl WH IM H?)0| : 
the contact time varies from the short collision time in 
the dynamic regime to th e much longer shear time in the 
quasi-static regime |l0lj . while the distribution of the 
force intensity is very wide. We shall not discuss in the 
following those distributions of contact time and force 
intensities. We shall rather focus on the following three 
quantities : coordination number, mobilization of friction 
and anisotropy of the contact forces. 



A. Coordination number 

As was shown on Fig. ^ the contact network is very 
sensitive to the inertial number. A small dilation of the 
material (around 10%) is enough to observe a transition 
from a dense contact network to multiple, or even binary, 
collisions between grains. A quantitative indicator is the 
coordination number Z, that is to say the average number 
of contacts per grain. The variations of Z as a function 
of / are shown in Fig.^J Z increases as I decreases, and 
tends toward a maximum value Z max when / — * 0. A 
possible fit is: 



(27) 
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which is drawn on Fig. which gathers the results for 
a given rigidity number n and various e and /i. 
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FIG. 13: Variation of the maximum coordination number as 
a function of k (e = 0.1, fi = OA). 




FIG. 12: Variation of the coordination number Z as a func- 
tion of I (k = 10 4 , e = 0.1 (*), e = 0.9 (o)) (a) ft = - 
Fits Z = 3.75 - 7.11/ ' 51 ( — Z = 3.00 - 6.04/ ' 41 (-) (b) 
li = 0.4- Fits Z = 2.84- 3.80/ ' 48 ( — ), Z = 2.60-5.48/ - 50 

(-)■ 



B. Mobilization of friction 

Inside the population of contacts, and for frictional 
grains, we now introduce a distinction between the "slid- 
ing" contacts where the local friction is completely mo- 
bilized (|T| = fj,N) and the other "rolling" contacts 
(\T\ < £tiV). This distinction is slightly different from the 
one proposed in |4jj > which distinguishes "fluid" contacts 
(collisions and sliding enduring contacts ) and "solid" 
contacts (rolling enduring contacts). We defin e Z R as the 
average number of sliding contacts per grain 103]. The 
Fig. E| (a) shows the variations of Z s as a function of the 
inertial number I. We observe that Z s increases with I 
in the quasi-static regime, up to a maximum in the inter- 
mediate regime. Moreover, the Fig.^](a) indicates that 
the Z S (I) curve depends on the restitution coefficient e. 

We have shown on Fig. ^] (b) the variations with / of 
the ratio M = Z s /Z, which, as the proportion of sliding 
contacts, is an indicator of the mobilization of friction. 
We observe that, contrarily to Z s , M increases, approxi- 
mately logarithmically, as a function of /. We also notice 
a slight increase of M when k increases. 



Fig. 1121 shows that the coordination number does not 
depend only on the geometry, through the solid fraction, 
but also on the mechanical properties of the grains e 
and \x. The exponent 7 is nearly constant (7 ~ i), but 
Zmax and c depend on e and ll. When jj, decreas es, Z 
increases and Z max tends to 4 for frictionless grains |lf)2j 
(see Fig. El (c)). We also notice that Z increases when 
e decreases, due to the increasing collision time. The 
Fig. E| indicates that Z max decreases significantly with 
Kj as expected. Among the various quantities which we 
have studied, the coordination number is the only one 
which varies significantly with the rigidity of the grains. 



C. Anisotropy of the contact network 

We now discuss the anisotropies of the contact net- 
work. We call (j> the direction of a contact counted 
counterclockwise from the x direction, between and tt. 
(n^^t,),) is the local frame in the direction <fi. Let us call 
(N) the average normal force in the homogeneous layer, 
then {N(4>)} and (T(<f>)) the average normal and tangen- 
tial forces in the homogeneous layer in the direction <j>. 
The anisotropies are described by the three angular dis- 
tributions of contact orientations p(4>), of intensities of 
normal forces £n(4>) — (N((f>))/(N) and of intensities of 
tangential forces £t(</>) = (T(4>))/(N). Then we define 
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FIG. 14: (Color online) (a) Variation of Z s as a function 
of I (fx = 0.4 - k = 4000 - e = 0.1 (M), e = 0.9 (o)) (b) 
Variation of M as a function of I, for various k,. 



CnW = p(4>)tiN(cb) and Ct(0) = p(4>)fr (</>)• Those an- 
gular distributions satisfy the normalization relations : 




0.02 
0.03 
0.04 




So p{4>)d(l) = l, 

< £ bfitW = 1, (28) 

k J^CtW^ = o. 

We show in Fig.Elthe two quantities CaK 1 ^) and Ct(^) 
which will be useful in the discussion of the friction law. 
We distinguish first quasi-static and dynamic regimes, 
second frictional and frictionless grains. A positive value 
of Ct(</>) indicates that the tangential forces induce on 
average a counterclockwise rotation of the grains, and is 
represented with white symbols in Fig. ^] (c). A neg- 
ative value of Ct(</>) indicates that the tangential forces 
induce on average a clockwise rotation of the grains, and 
is represented with black symbols in Fig. 1151 (c) . 



FIG. 15: (Color online) Angular distribution of the contact 
forces (quasi-static regime (I — 0.005 - O), dynamic regime 
(I = 0.13 - o): (a) 0v(<« (p = OA), (b) Ov(0) (p = 0), (c) 
Ct(4>) (p. = OA) - > white symbols, < black symbols. 



We notice a strong anisotropy of the contact network 
with privileged orientations for Cn{4>) along the direc- 
tions of shear (0 ~ and 7r) and of maximum compres- 
sion (<j) ~ 27r/3), and for Ct(</>) along the directions of 
shear (</> ~ and 7r) and of the shear gradient (cf> ~ 7r/2). 
Those anisotropies slightly change between the quasi- 
static and dynamic regimes and between frictional and 
frictionless grains cases. 

Those anisotropies may be explained within a very sim- 
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plified picture of a granular material organized in layers 
along the shear direction. Then, there are two kinds of 
contacts between grains, inside a layer (<fi ~ and it) and 
between layers (tt/3 < </> < 2tt/3). The contacts between 
layers are created along the direction of maximum com- 
pression ((f) ~ 2tt/3). In the quasi-static regime, those 
contacts are maintained up to the point where the grains 
separate (<f> ~ tt/3). This is in contrast with the dynamic 
regime where the grains bounce, so that is stronger 
around 7r/3 and smaller between 7r/3 and 27r/3. 

We observe (see Fig. ^] (a)) that the tangential 
anisotropy is well described by the following expression : 

CtW = /t(J)Cw(0)cos(2(«, (29) 

with an increasing positive function /t(-0- This means 
that the contacts between layers favor a clockwise rota- 
tion of the grains, while the contacts inside layers favor 
a counterclockwise rotation of the grains. Furthermore, 
this shows that the average tangential force is smaller in 
the quasi-static regime than in the dynamic regime. 

VIII. MICROSCOPIC ORIGIN OF THE 
FRICTION LAW 

We now try to understand quantitatively the macro- 
scopic friction law on the basis of the previous micro- 
scopic information (fluctuations and anisotropy of the 
contact network). 

A. Friction and rotation 

We first discuss the tangential anisotropy and show 
how it is related to the rotation of the grains 98] . 

1. Average rotation 

In a first step, we take into account the average rota- 
tion velocity cu. We consider the homogeneous shear of an 
assembly of grains of diameter d with an average shear 
rate j. Then, the tangential relative velocity between 
two grains (see Sec. Ill Bp is given by : 

V T ((f>) =d(u + j sin 2 ((/))). (30) 

The average tangential force exerted on a grain should 
be equal to zero in steady state : 




pwntw = o. (3i) 



We start with two very crude assumptions : all the con- 
tacts are sliding (T = —fj,Nsign(V T )) and the normal 




I 



FIG. 16: Angular distribution of tangential contact forces: 
(a) Ct(<P)/Cn{4>) (°) fitted by f T cos(2^) (-) for I = 0.13 
and /i = 0.4, (b) f T (I)/n for fi = 0.4 ■ ) and fi = 0.8 (o) 
compared with ^-I 1 ^ 3 . 



force distribution is isotropic (Ov(0) = 1/tt)- Conse- 
quently : 



/ sign(u + j sin 2 (<j)))d(j) = 0. (32) 
Jo 

This is an explanation for the relation ui = —^7, which 
was described in Sec. IIVB II (see Fig. |21 (b)). Con- 
sequently V T ((f)) = —d A /cos(2(f>)/2, Qt — m/tt when 
< 4> < 7r/4 and 37r/4 < (f> < tt and £t = ~n/ 7r when 
7r/4 < (j> < 37r/4. In comparison with Fig. El( c )> we no- 
tice that the sign of Cr(</0 is correct, but that the order 
of magnitude is too large, by a factor of around 10. 



2. Fluctuations of the rotation 

As a way to understand the order of magnitude of 
C/r(<^), we now take into account the fluctuations of the 



19 



rotation velocity, which have been evidenced in Sec. IVII 
Denoting as 5u>i j the fluctuations of rotation of two 
grains i and j in contact, their relative tangential ve- 
locity at the contact point becomes : 

V£(<j>) = -d7Cos(20)/2 + d/2{6uj l + Sluj). (33) 
Keeping the assumption of sliding contacts, we predict : 



quasi-static regime, we have observed that most of the 
contacts are rolling. Furthermore, our simulations reveal 
correlations of the rotations of grains in contact: it seems 
that the flowing granular material is organized in clusters 
of grains rotating in the same way |98j . Such correlations 
of the grain motion deserves further study [89l l95| . 

B. Friction and anisotropy 



Ct(</>) =-K N ((b)(sign(V T )}, 



(34) 



where (sign(V T )} is the statistical average over the fluc- 
tuating rotations of the two grains. Discrete numerical 
simulations have shown that the distribution of the rota- 
tion velocity is approximately lorentzian |55| . We make 
the assumption that the fluctuations of rotation of two 
grains in contact are not correlated. Then the random 
variable (Suji + Sujj)/2 follows a lorentzian distribution, 
with a zero mean value and a variance Scj/V2. Then we 
obtain : 



(sign(V T )) = arctan \ —X — cos(20) I . (35) 

7T \V2Su; 

Using the rotation scaling law (|24(1 : 



2 /I 1 / 3 
(sign(V T )} — arctan ( — y=- cos(20) ) , (•'>(>) 



and for small / (< 0.2) : 



We now discuss the friction law on the basis of 

the information on the contact network. We would like to 
understand the increase of /i* with /, and its dependence 
with the microscopic friction /i. 

A first possible interpretation lies in the increase of 
the mobilization of friction M(J) (Fig. [21(b)): most of 
the contacts are rolling in the quasi-static regime, while 
most of them are inelastic sliding collisions in the dy- 
namic regime. Consequently, the energy dissipation, and 
hence the effective friction, should be stronger in the dy- 
namic regime. However, since the effective friction coef- 
ficient of an assembly of frictionless grains is not equal 
to zero, this interpretation is certainly not sufficient. We 
are now going to show the crucial role of the anisotropics 
of the contact network. 

We consider the homogeneous shear of an assembly of 
grains of diameter d with average solid fraction i>, coor- 
dination number Z, and normal force (N). The stress 
tensor is dominated by the contribution of the contacts 
(see Eqn. I12|) . It is possible to express it as a function 
of the angular distributions of c ontact fo rces, which we 
have previously defined [O, EE llCHl 11051 ] : 



(signiV 1 )) 



T ))^-— I 1/3 cos(2^), 



so that at the end : 



V2 



CtW^ — m/ 1/3 6vWcos(20). 



(37) 



(38) 



This expression reproduces the observed angular de- 
pendence in Cn{4>) cos (2(f)), shown in Fig. ^| (a). The 
prediction for the dependence of the amplitude on / and 
/i is in agreement with the observed pre- factor /t(-0 de- 
fined in Eqn. (|29() . as shown in Fig. Q2l(b). The order of 
magnitude is now consistent with the observations. 

As a conclusion, the fluctuations of the rotation veloc- 
ity are a possible quantitative explanation of the angu- 
lar distribution of tangential forces Or (</>)• When / de- 
creases, the relative fluctuations of rotation increase, so 
that the average relative tangential velocity of two grains 
in contact tends to zero, which kills the frictional effect. 
This model is very crude. In the dynamic collisional 
regime, we have seen that most of the contacts are sliding 
(see Fig. I14f) and the assumption of uncorrelated fluctu- 
ations may seem reasonable. On the contrary, in the 



S = 



2vZ < N > 



nd 



[CaK0)™0 + Ct(0)^] «> ^d<?f>. (39) 



Using the properties of the stress tensor (Yi xy — 
and Y, xx = "S yy ) and the normalization of angular distri- 
butions, the effective friction coefficient takes the simple 
expression : 



[6v($ sin(20) + Ct W cos(20)]#. (40) 



In the first term associated to the normal forces (fi%), 
the factor sin(2</>) is positive for (j> between and tt/2, 
so that it decreases the effective friction, and negative 
for (j> between tt/2 and ir, so that it increases the effec- 
tive friction. Consequently, the evolution of the angular 
distribution between the quasi-static regime and the dy- 
namic regime (see Fig. llSl fa) and (b)) might explain part 
of the increase of the effective friction. 

It is possible to give an estimation of the second term 
associated to the tangential forces (n^), using the ap- 
proximation for Ct(^), with CaK0) = I/ 71 "- We ob- 
tain /Xy ~ — (/i/v^tt)^ 1 ^ 3 - This contribution is negative 
and small. The complete calculation (Fig. El (b)) con- 
firms that the contribution of tangential forces to the 
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effective friction remains of the order of 15% 

for n = 0.4 and 30% for /i = 0.8. 

Fig- 133 (a) compares the prediction based on the angu- 
lar distribution l|40l) with the complete calculation based 
on (|12|l . both for frictional and frictionless grains (the 
results of Fig. UJ). The agreement is very good for fric- 
tionless grains, but we observe a deviation for frictional 
grains, of the order of 15%, which we do not explain at 
this stage. A possible explanation is the small polydis- 
persity since the monodispersity was the only simplifying 
assumptions leading to (|4U|) . 

This result is very paradoxical. Contrarily to our first 
tentative explanation of the friction law, rather than ex- 
plaining the macroscopic friction, the microscopic friction 
(tangential forces) has a small and even negative con- 
tribution to the macroscopic friction. The friction law 
depends mostly on the angular distribution of normal 
forces. A small variation of the distribution is enough 
to increase the effective friction by a factor of two be- 
tween the quasi-static and the dynamic regime. In the 
same way, this distribution is more isotropic in the case 
of frictionless grains, and this explains the shift of the 
friction law. The microscopic friction has an indirect ef- 
fect on the friction law, through the modification of the 
angular distribution of normal forces. This would mean 
that the very origin of the visco-plastic constitutive law 
relies in the anisotropy of the contact network (fabric) 
in response to the shear. This point which has already 
been studied in the quasi-static regime |13| would deserve 
further study. 
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IX. INFLUENCE OF GRAVITY 

Up to now, we have studied the simplest flow geome- 
try, plane shear between two rough walls, without gravity. 
Then the stress distribution is homogeneous. In the limit 
of rigid grains and in the intermediate regime, we have 
observed steady homogeneous flows, in which we have 
measured the dilatancy and friction laws. Other simu- 
lations of homogeneous shear stat es without walls have 
confirmed those observations |l06j| . 

We now wonder if those laws are valid in other flow ge- 
ometries (annular shear, vertical chute, inclined plane...) 
where, as a consequence of the heterogeneous stress dis- 
tribution, the flow is no more homogeneous, and often 
localized near the walls or the free surface |25l l2q | . 

As a simple example of an heterogeneous stress dis- 
tribution, we still consider the plane shear flow geome- 
try, but we study the influence of gravity g along y (ori- 
ented toward the fixed wall). Then, the shear stress S 
remains constant in the sheared layer, but the pressure 
varies along y according to : 



P(y) = P w + Pg g f v{t)dt, 
Jy 



(41) 



FIG. 17: (Color online) Friction and anisotropy : (a) Com- 
parison of Eqn. (black symbols) with direct measurement of 
Fifl-Q (°P en symbols) (fi = (D), n 7^ (o)), (b) Tangential 
contribution for /j, — 0.4 ■ ) and (1 = 0.8 (o ). 



We notice that this kind of flow has b een p reviously 
studied experimentally |63l l64l l65l l83l Il07| numeri- 
cally [H Hfl IH E3 and theoretically |23 . 



where P w is the pressure prescribed on the moving wall. 



A. Dimensional analysis 



We define dimensionless numbers from the various time 
scales of the system, and notice that in an heterogeneous 
system, those quantities depend on y. In addition to the 
shear time 1/7(2/), the inertial time yj m/P(y), and the 
collision time r c , the presence of gravity introduces an- 
other gravity timescale yj d/g. Consequently, in addition 
to the inertial I(y) and rigidity n(y) numbers, we must 
consider a third dimensionless number, measuring the 
intensity of gravity. We call G(y) this gravity number, 
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defined as the ratio of inertial to gravity times : 



G(y) 



9™ 

P(v)d' 



(42) 



In the following, we have not studied the influence of 
the rigidity of the grains : n(y) varies between 10 4 near 
the moving wall and 10 2 near the fixed wall, where we 
are beyond the rigid grain limit. The simulated systems 
{H/d ~ 40) are described by the two dimensionless num- 
bers I g = V/H \JmjP w and G = J -jr^- Ig varies be- 
tween 0.0025 and 0.05, whereas G varies between and 
3.2. 




B. Localization 

Let us first assume that the dilatancy and friction law 
identified in homogeneous shear flows are still valid in 
presence of gravity. Then we predict that, like S/P{y), 
I(y) should decrease from y = H to y = 0. Consequently 
the shear rate should decrease and the solid fraction in- 
crease. For a large enough G, there is a point (H — A) 
where S/P(H — A) = [i* min - Below this point, the granu- 
lar material should be at rest: 7 = and v — v max . This 
would mean that the shear is localized in a layer of width 
A. 

The measured velocity profiles v x (y) (Fig. ED ( a )) in- 
dicate a shear localization near the moving wall. Those 
profiles are approximately exponential in the upper part 
(see v x(y) — ^ ex P(— -^x^)i Ir0m which we esti- 

mate the width of the shear layer A. Contrarily to the 
annular shear or vertical chute geometries [25| , where this 
length is always of the order of five to ten grain diame- 
ters, this length strongly varies as a function of G (Fig. 1181 
(b)) and diverges when G — > 0, which means that the flow 
becomes homogeneous in the limit of small gravity. The 
solid fraction increase at distance from the moving wall 
and is very high in the quasi-static region. 
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FIG. 18: (Color online) Influence ofG: (a) Velocity profiles 
(logarithmic scale) (b) Width X of the shear layer (H/d = 
40, e = 0.1, fi = 0.4, I g = 0.015; 



translation and rotation scaling laws l|24|l remain valid 
but the exponents a and (3 increase respectively from 0.5 
to 0.6 and from 0.3 to 0.5. 



C. Dilatancy and friction laws 

For each value of G, we deduce from v(y), 7(2/), P{y) 
and S(y) the variations of the local solid fraction and of 
the effective friction coefficient as a function of the local 
inertial number (Fig. ^5] (a) e t (b)). A single simulation 
provides the whole curve, but the data are noisy since 
they result from local measurements, which are not av- 
eraged over the whole layer. The tendency is analogous 
to the no-gravity case: decrease of v from v max , and in- 
crease of /x* from \x* min when / increases. However, as 
G increases, v max increases while \x* min decreases. The a 
and b factors are affected as well. 

We have also measured the influence of G on the con- 
tact network and on the fluctuations [2(| : G has nearly 
no influence on the Z(I) and M(I) dependencies; the 



D. Comments 

Those measurements indicate that the measurements 
performed in homogeneous shear flows are not sufficient 
to describe other flow geometries. Then, the inertial 
number is not the only relevant quantity to describe the 
shear flow. We must take into account other quantities, 
such as the gravity number in plane shear flows with 
gravity. At this stage, we are not able to explain quan- 
titatively the observations which we have briefly quoted, 
such as the nearly exponential shear localization, or the 
decrease of \i* m m with G. Following our explanation of 
the friction law in the homogeneous case (see Sec. IVIli)l . 
we suggest that this could be related to the evolution of 
the anisotropy of the contact network: as G increases, 
the distribution (n(4>) would increase for < (f> < ir/2. 
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FIG. 19: (Color online) Influence ofG: (a) Dilatancy law (b) 
Friction law (logarithmic scale) (e — 0.1, n — 0.4, I g = 0.015 J. 



X. CONCLUSION 



values of / correspond to the quasi-static critical state 
regime of soil mechanics, while large values of I corre- 
spond to the dynamic collisional regime of the kinetic 
theory. When / increases in the intermediate regime, 
we have measured an approximately linear decrease of 
the solid fraction from the maximum packing value, and 
an approximately linear increase of the effective friction 
coefficient from the static internal friction value. From 
those dilatancy and friction laws, we have deduced a 
visco-plastic constitutive law, with a plastic Coulomb 
term and a viscous Bagnold term. We have also mea- 
sured scaling laws for the relative velocity fluctuations 
as a function of /. We have shown that the mechani- 
cal characteristics of the grains (restitution, friction and 
elasticity) have a very small influence in the intermediate 
regime. We have shown that the origin of the friction law 
relies in the fluctuating rotations of the grains and on the 
angular distribution of contact forces. 

The study of plane shear flows with gravity has 
shown that the generalization to other heterogeneous 
flow geometries is not straightforward. However, the 
study of annular shear flows j2|J and inclined plane 
flows [M El E3 reveal the 

same qualitative tenden- 
cies. Other flow geometries such as vertical chute, ro- 
tating drum and heap flow should be analyzed in the 
same way 25;] . 

We notice that a generalization of those ideas to steady 
uniform shear flows of cohesive granular materials has 
been successful |106| . 

In this paper, we have restricted our attention to ve- 
locity controlled shear flows, so that it was not possible 
to study the flow threshold. We think that it was ev- 
idenced indirectly through the appearance of intermit- 
tencies for small enough / j9^. A specific study of the 
jamming mechanisms should be performed b y co ntrolling 
the shear stress, either in plane shear flows [4(| or down 
inclined planes [2a. r71|. 



We now summarize our conclusions. We have consid- 
ered the simplest flow geometry (plane shear without 
gravity), where the stress distribution is homogeneous. 
Using molecular dynamics simulation, we have submit- 
ted a dense assembly of frictional, inelastic disks, to a 
given pressure and shear rate. We have observed steady 
homogeneous shear flows, which become intermittent in 
the quasi-static regime. We have shown that, in the limit 
of rigid grains, the shear state is determined by a single 
dimensionless number, called inertial number /, which 
describes the ratio of inertial to pressure forces. Small 
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